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Abstract 

The process is calculated to two- loop order (p 6 ) in Chiral Perturbation Theory 
(ChPT) in the isospin conserved case. We present expressions suitable for use with 
previous work in two-loop CHPT where the order p 4 parameters {L\) were determined 
from experiment. We point out that all the order p e parameters (Cf ) that appear 
in the value of /+(0) relevant for the determination of \V US \ can be determined from 
K13 measurements via the slope and the curvature of the scalar form-factor. As by 
product we update the value of the CKM matrix element | V us \ . 
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1 Introduction 



Weak semileptonic kaon decays to a pion and a lepton-neutrino pair {K^) have a long 
history. The early theoretical treatments can be found in the review [1] . This decay plays an 
important role in the determination of the CKM matrix-element V us , see e.g. the discussion 
in [2] or [3]. The theoretical basis for this determination is the paper by Leutwyler and 
Roos [4]. The basis for this evaluation were radiative corrections and estimates of one- loop 
chiral corrections. The full Chiral Perturbation Theory calculation to order p 4 (see Sect. 2 
for a short explanation) was performed by Gasser and Leutwyler [5]. References to earlier 
work on the non-analytic corrections can be found there. Recent Reviews of the situation 
can be found in [6] or [7]. 

Since that time a lot of work has been performed in Chiral Perturbation Theory. An 
update of the calculation of [5] to order p 6 is thus necessary. Partial studies are done, 
the double logarithm contribution is small [8] and a possibly large role for terms with two 
powers of quark masses has been argued for in Ref. [9]. One full order p 6 calculation exists 
[10], but it uses outdated values of the ChPT constants as well as an older version of the 
classification of p 6 constants. In this paper we present an independent calculation of the 
Kij, amplitudes to order p 6 in ChPT in the isospin limit. We present numerical results 
with values for the ChPT constants resulting from fits to order p 6 [11, 12, 13, 14]. Related 
work is the update of the electromagnetic radiative corrections given in [15]. 

We present a few definitions of ChPT in Sect. 2 to determine our notation. Sect. 3 
defines the form-factors used in Kg 3 decays. Analytical results are presented explicitly 
in Section 4 for the form-factors up to order p 4 and for the part depending on the order 
p 6 parameters {CI). The remaining parts are rather long and can be obtained from the 
authors on request. Section 5 describes one of our main results, the fact that all needed p 6 
constants for the value of /+(0) can be experimentally determined from experiments. 
The value of /+(0) is of course needed for determinations of V us and is of use for future 
precise measurements of K — > -nvv. 

In Section 6 the presently available data set is discussed. Here we also point out that 
the often used linear approximation in the form-factors can have a sizable effect on the 
measured value of the slope and the value at t — 0. This effect is of similar size as the 
experimental errors quoted. We present an extended discussion of the numerical results in 
Sect. 8 after a short discussion of the inputs used in Sect. 7. Our final conclusions for A + 
are in Sect. 8.3 and of /+(0) in Sect. 8.5. We summarize our results in Section 9. 

2 Some definitions 

Chiral Perturbation Theory is the modern way to derive the predictions following from the 
fact the SU(nf)i x SU (^/)_r chiral symmetry in the limit of rif massless flavours in QCD is 
spontaneously broken by nonperturbative QCD dynamics to the diagonal vector subgroup 
SU{rif)v- It is the effective field theory method to use this property at low energies. 
It takes into account the singularities associated with the Goldstone Boson degrees of 
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freedom caused by the spontaneous breakdown of chiral symmetry and parametrizes all the 
remaining freedom allowed by the chiral Ward identities in low energy constants (LECs). 
The LECs are the freedom in the parts of the amplitudes that depend analytically on the 
masses and momenta. The expansion is ordered in terms of momenta, quark masses and 
external fields. Recent lectures introducing this area are given in ref. [16]. We use here the 
standard ChPT counting where the quark mass, scalar and pseudoscalar external fields are 
counted as two powers of momenta. Vector and axial- vector external currents count as one 
power of momentum. The lagrangian can be ordered as 



••effective 



£ 2 + £ 4 + £ 6 + ■ ■ ■ 

10 90 94 



i=91 



The index i in £j stands for the chiral power. The precise form of £ 2 and £4 is given below 
while £ 6 can be found in [17]. The lowest order, 0(p 2 ), in the expansion corresponds to 
tree level diagrams with vertices from £ 2 , the next-to-leading order, NLO or 0(p 4 ), to 
one-loop diagrams with vertices from £ 2 or tree level diagrams with one vertex from £ 4 
and the rest from £ 2 . The next-to-next-to-leading order, NNLO or 0(p 6 ), has two-loop 
diagrams, one-loop diagrams with one vertex from £ 4 and tree level diagrams with one 
vertex from £ 6 or two vertices from £4 and all other vertices from £ 2 . The loop diagrams 
take all singularities due to the Goldstone Bosons correctly into account. The singularities 
are the real predictions of ChPT while the other effects from QCD are in the values of 
the LECs. The diagrams, in addition to wave-function-renormalization, relevant for the 
processes discussed in this paper are shown in Figs. 1, 2 and 3. 

The expressions for the first two terms in the expansion of the Lagrangian are given by 
(F refers to the pseudoscalar decay constant in the chiral limit) 

£2 = ^-{(D^D^U) + ( X t[/ + x [/t )} > (2 .2) 



and 



£4 = U (D^D^U) 2 + L 2 {DprfDJJ) (D^D U U) 

+L 3 (D^D fl UD u U^D u U) + U(D^D^U){x ] U + X U ] ) 
+L 5 (D^D tl U(x i U + rf X )) + L 6 (x^U + X ^) 2 
+L 7 ( X j U - x^f + L S { X WU + X^xrf) 

-iL 9 (F*D ll UD u U t + Fj; v D»rfD v U) + L W {U ] 'F*U F L ^) , (2.3) 



while the next-to-next-to-leading order is a rather cumbersome expression [17]. The special 
unitary matrix U contains the Goldstone boson fields 



U = exp Cj^m\ , M = 



( >° + 7^ - + K+ \ 
K- K° 



(2.4) 



2 



The formalism is the external field method of [18] with s, p, 1^ and matrix valued scalar, 
pseudo-scalar, left-handed and right-handed vector external fields respectively. These show 
up in 

X = 2B (s + ip) , (2.5) 

in the covariant derivative 

DpU = d„U - irJU + iUlp , (2.6) 

and in the field strength tensor 



F^=dAr)u-d v l(r) lt -i[l(r) lt ,l(r) v ] 



For our purpose it is sufficient to set 



m 



m 



m s j 



I 

' y/2 



V U(j W+ V W W+ \ 



v: d w- 
\ v: s w- 



with g 2 the weak coupling constant, related to the Fermi constant by 



9l 



V2 Smlr 



(2.7) 



(2.8) 



(2.9) 



2.1 Renormalization Scheme 

We use the renormalization scheme as explained in [17] and [19]. It extends the scheme 
from [18] naturally to two-loops. Notice that the work of Post and Schilcher [10, 20, 21] 
used a slightly different scheme. The scheme employed here does not introduce the e 2 term 
in Eq. (39) of Ref. [10]. Subtractions are performed via 



with 



Ci = (^) 2(d - 4) (q(/^)-^[rp ) A 2 + (rf) + rf)( / ,))A 



lnC = 4(ln(4.) + r'(l) + l), A = T ^I) 



(2.10) 



(2.11) 



The coefficients I\ can be found in [18, 17] and the vf \ T^rfV) in [17]. We will in the 
remainder always write CI and L\ but the \i dependence is of course present. 



3 The Kis form- factors: definition and 0(p 4 ) 

The decays we consider are 



K + {p) 
K°(p) 



■K~{p')t{pi)Vt{Vv) , 



(3.1) 
(3.2) 
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and their charge conjugate modes, i stands for /i or e. 

The matrix-element for K^, neglecting scalar and tensor contributions, has the struc- 
ture 

T = ^V:^F, + (p',p), (3.3) 

with 

F, + (P,P) = <n°(p')\V^(0)\K + (p)>, 
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= ^[(?/+pU+ (*) + (p-A/- (*)]■ (3-4) 
To obtain the K® 3 matrix-element, one replaces by 

F»°(p',p) = <n-(p>)\V^ 5 (0)\K°(p)> 

= (p'+p)jf^(t) + (p-p'),f K ^-(t). (3.5) 

The processes (3.1) and (3.2) thus involve the four K m form-factors /f +7T °(t), f£° n (t) 
which depend on 

t=(p'-p) 2 = (pi+Pu)\ (3.6) 

the square of the four momentum transfer to the leptons. 
In this paper we work in the isospin limit thus 

J± = J± = J± = J± ■ (3.7J 

jKn j g re f errec i to as the vector form-factor, because it specifies the P-wave projection of 
the crossed channel matrix-elements < | V^~ l5 (0) | K + ,7i° in >. The S-wave projection 
is described by the scalar form-factor 

fo(t)=f + (t)+ t /-(*)■ (3.8) 
m K — 



Analyses of data frequently assume a linear dependence 



U,o(t) = / + (0) 



t 

1 + A+,o J" 

ml.. 



(3.9) 



For a discussion of the validity of this approximation see [5] and references cited therein. 
We will discuss it to order p 6 and in comparison with the data. At the expected future 
precision it will be necessary to go beyond this approximation. 
Eqs. (3.8) and (3.9) leads to a constant f-(t) , 

/-(*) = /-(0) = / + (0)(A - \ + ) ml ~ ml . (3.10) 

m 7T+ 
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(a) (b) (c) 

Figure 1: The tree level diagrams contributing to the form- factor. A • indicates a 
vertex from C 2 , a X a vertex from £ 4 and a (g> a vertex from £ 6 . The insertion of the weak 
current is indicated by the wiggly line. 

The form-factors f±,o(t) are analytic functions in the complex t-plane cut along the positive 
real axis. The cut starts at t = (m.# + m 7r ) 2 . In our phase convention, the form-factors are 
real in the physical region 

nrf <t< {m K — m n ) 2 . (3-H) 
A discussion of the kinematics in K$% decays can be found in [6] . 

4 Analytical Results 

The total result we obtain is split by chiral order. 

m) = f?\t) + fl A \t) + ff\t), (« = +,- o). (4.i) 

4.1 Order p 2 

This has been known for a very long time and is fully determined by gauge invariance. 

/£ ) (t) = /o (2) (f) = l,, f ( -\t) = 0. (4.2) 

It results from the diagram in Fig. 1(a). 

4.2 Order p 4 

This was first calculated within the ChPT framework by Gasser and Leutwyler in 1985 [5] . 
The result contains the nonanalytic dependence in the symmetry parameters predicted by 
[22]. The form of the form-factors we use (which is equivalent to the result of [5] to order 
p 4 ) is the one which our expressions for the p 6 contribution correspond to. 



Flf^it) = 2L r 9 t + 3/8 (A{m 2 v ) + A{ml) + 2A(m 2 K )'^ 

-3/2 (B 22 (ml,m 2 K ,t) + B 22 (m 2 K ,m 2 v ,t)) . (4.3) 
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F 2 f^(t) = (4m 2 KZ 4ml)Ll-2(ml-ml)Ll + l/2A(ml)-5/12A(ml) 
+7/12A(m 2 K ) + B{ml, m 2 K , *)(-l/12m^ - 5/12m 2 K + 5/12t) 
+fi(m|,m;,t)(l/12m2 - 7/12m^ + 1/4*) 
+Si(rr4, m^, t)(-7/12rr4 + 19/12m^ - 5/12*) 
+B 1 (m 2 K ,m 2 ,t) (-11 /12ml + 23/12m^ - 1/4*) 
+5 21 (m2,m^,t)(3/2m2 - 3/2m 2 K - 5/6*) 
+5 21 (m^,m2,t)(3/2m2 - 3/2m 2 K - 1/2*) 

+5 22 (m2 , m^, *)(-5/6) + 5 22 (™L mj, *)(-l/2) . (4.4) 

These results are obtained from the diagrams in Figs. 1(b), 2(a) and 2(c), together with 
wave-function renormalization. 



4.3 Order p 6 

The p 6 contribution we split in several parts 

tf\t) = j- 4 (fC(t) + f t L (t) + f t B (t) + f t H (t) + fr(t)) , (i = +,-,0). (4.5) 

7T 

The split between the last three terms is not unique and depends on how the irreducible 
two-loop integrals are separated from the reducible ones. The first two terms are the ones 
containing the dependence on the p 6 and p 4 coupling constants. 

The ones with dependence on C[ stem from wave-function renormalization and the 
diagram of Fig. 1(c). The results are 



fi(t) = R*Z + tR«? + t 2 R«Z 

, 2 



-8 (ml - 



K 



(C[ 2 + CJJ + t - 4ml(2C{ 2 + 4C[ 3 + C r M 



m 



-8G 



12 



32C[ 3 



8^63 



+ t 2 (-ACl, + AG 



and 



90; > 



+2 r>Kn 



(4.6) 



/_(*) = R-o + tR_i + 1 R 

= (m\ - ml) [ml (+24C[ 2 - 16C[ 3 + 8C[ 5 + 16C[ 7 + 8C 3 r 4 + 4C t 



+4Q 5 + 4CS 



90 ; 



m 



24CT, + 32C[ 3 + 16C[ 4 + 16CT R + 8C : 



+AC r m + 8C fi r 4 + 4C, 



64 



"90V 



'12 

*4 (m^ 



14 



'15 



7 34 



2C[ 2 + a r 



90 J 



(4.7) 



We have followed a notation very close to the one in [14]. Notice that we have the 
relations 



R K _l = 0, 
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(a) (b) (c) (d) 



(e) (f) (g) (h) 




Figure 2: The one and two- loop diagrams without overlapping loop integrals contributing 
to the Kf^ form- factor. A • indicates a vertex from £ 2 > a X a vertex from £ 4 and a <S> a 
vertex from £ 6 . The insertion of the weak current is indicated by the wiggly line. 
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i 




(b) 



Figure 3: The two-loop diagrams with overlapping loop integrals contributing to the 
form-factor. A • indicates a vertex from C 2 , a X a vertex from £ 4 and a a vertex from 
Cq. The insertion of the weak current is indicated by the wiggly line. 

z>Kw pi" f?^ + 



K +1 



\ (RU + + < . (4.8) 



The other are similar combinations of the C\ but in the expansion of the electromag- 
netic form-factors [14]. Notice that the last relation is really Sirlin's relation [23] and the 
second satisfies it as well. 

We have not quoted the remaining formulas, but will quote below some approximate 
numerical expressions. The exact formulas can be obtained from the authors on request. 
Our expressions satisfy the Ademollo-Gatto theorem [24]. 



5 Getting the value of /+(0) 



One of the problems we face here is whether the needed C\ can be determined from exper- 
iment. There are many of these coefficients showing up but as is obvious from Eq. (4.6), 
what we need is a value for C'[ 2 + £34- It turns out that this combination can actually 
be determined from measurements. The derivation given below relies on the fact that 
we need values for the p 4 constants, determined to order p 4 only, in the order p 6 part to 
be correct to the accuracy that we are working. We can determine all needed L\ to this 
accuracy relying only on data. 
We construct the quantity 

f (t) = f+(t) + 2 1 2 (/_(<) + 1 - F K /F n ) = f (t) + 2 t - (1 - F K /F V ) . (5.1) 

This has no dependence on the L\ at order p 4 , only via order p e contributions. Inspection 
of the dependence on the G\ shows that 

-^Cl 2 + A(t)+A(0). (5.2) 
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We emphasize that the quantities A(t) and A(0) can in principle be calculated to order 
p 6 accuracy with knowledge of the L\ to order p 4 accuracy. In practice, since a p 4 fit will 
include in the values of the L\ effects that come from the p 6 loops (due to the fitting to 
experimental values) we consider the p 6 fits to be the relevant ones to avoid double counting 
effects. Numerical results will be discussed in Section 8. 

The definition in (5.1) has essentially used the Dashen-Weinstein relation [25] to remove 
the LI dependence at order p 4 . It has also the effect that it removed many of the C\ from 
the scalar form-factor as well. The corrections which appear in the Dashen-Weinstein 
relation are include in the functions A(t) and A(0), these have both order p 4 [5, 22] and 
order p 6 contributions. 

It is obvious from Eq. (5.2) that the needed combination of C[ can be determined from 
the slope and the curvature of the scalar form-factor in decays. 

It seems possible that C[ 2 can be measured from the curvature of the pion scalar form- 
factor near [26]. When this calculation is complete, one can use the dispersive estimates 
of the pion scalar form-factor together with only a A measurement in K^ 3 to obtain the p 6 
value for /+(0). There are also some dispersive estimates for the relevant scalar form-factor. 
Unfortunately, these are not in a usable form at present [27]. 

The discussion of Ref. [9] can be shown in this light too. The constant A 3 of gener- 
alized perturbation theory correspond to a combination of the C\ from [17]. The precise 
combination is 



1 1 

"°34 "+" 2 v°14 + °17 "+" U 26 + °29 "+" °3lJ "+" g°33 



(5.3) 



As we have shown it is possible to eliminate all but C^ 4 , so this plays the role of A 3 here. 
C[ 2 is higher order in the counting employed in [9] and was not considered there. 

In Ref. [9] the additional observation was made that a relatively small change in the 
ratio Fk/F^, together with an adjustement of the constant A 3 can accommodate the CKM 
unitarity. 



6 Data 
6-1 

The most useful data points are those where the full dependence on the kinematical variable 
t = q 2 /m1 is shown. That means that the experiments that determined the value of Ao 
from the branching ratio or did not provide an actual t dependence but just fitted the 
linear slope in a global way are not that useful for us, but see below. 

Of the more recent experiments that quote data not from the branching ratio, the ones 
that gave a plot or numbers for f (t) are [28, 29, 30, 31]. Ref. [32] gives a plot but mentions 
that it is not statistically significant, which inspection of the plots confirms. 

There are some obvious problems with the data. E.g. the fo(t) from [28] do not go 
to 1 at t — 0. We have shown these data in Fig. 4 together with a linear and a pole 
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Figure 4: The data on /„(t)// + (0) from X^. The data are Clark77 [29], Donaldson74 [28], 
Buchanan75 [31] and Birulev81 [30], The latter reference has two distinct data sets. For 
comparison a linear approximation and a pole approximation with a mass of 800 MeV are 
shown as well. 



approximation corresponding to a mass of 800 MeV. This shows the accuracy needed to 
see the curvature. 

6-2 k+ 3 

Here we have not been able to find data that show a plot of fo(t). All experiments 
are analyzed in terms of a constant form-factor as discussed in [2]. There is one more 
experiment [33] that quotes measurements of Ao not included in [2]. 

6.3 K% 

Here the data are dominated by the recent high statistics CPLEAR data [34]. There 
exists a very high statistics older experiment [35]. They provide plots with different data 
assumptions and can thus not be easily compared at the level of f+(t) directly. But [35] 
quoted both a linear and quadratic fit to f+(t). In order to show the relation between the 
most recent and older data we have plotted the data of [31], [30] and [34] in Fig. 5. 
We have performed some simple fits of the form 

/+(*) =a+(l + X + -^- + c + A (6.1) 
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form 


a + 


A+ 


c+ 
[GeV" 4 ] 


10 5 

[GeV" 2 ] 


10 3 


En (6 1 "l 
Eq. (6.1) 
Eq. (6.1) 
Eq. (6.1) 


= i 
1.000 ±0.004 

= 1 
1.008 ±0.009 


0245 + 0006 
0.0245 ±0.0015 
0.0238 ±0.0017 
0.0181 ±0.0068 


= 

= 
0.5 ± 1.2 
2.8 ±2.8 






Eq. (8.2) 
Eq. (8.2) 
Eq. (8.2) 
Eq. (8.2) 


1.008 ±0.008 

= 1 

1.009 ±0.004 


0.0180 ±0.0067 
0.0236 ±0.0019 
0.0201 ±0.0006 
0.0170 ±0.0015 


2.7 ±3.0 
0.4 ± 1.2 

3.2 

3.2 


-4.3 ±2.5 
-2.2 ±0.7 
-4.7 ±0.5 
-4.7 ±0.5 


0.19 ±0.21 
0.02 ±0.09 
= 0.22 
= 0.22 



Table 1: The fits to the CPLEAR data of various theoretical forms of / + . In the last four 
fits, which use th ChPT results, A + and c + are derived quantities. The symbol = means 
this quantity was set to this value in the fit. 



to the CPLEAR data. The fits agree extremely well with those reported in [34] and are 
given in Table 1. Notice that the fits that go beyond the linear approximation and leave 
the normalization free give a significantly lower A + and with larger errors. It is within its 
errors compatible with the linear fit, but outside the errors from the linear fit, we consider 
that result from the fit with curvature to be more reliable. The shift is of similar size to 
that observed in Table 1 in [35]. 



6-4 K+ 

There are recent high statistics experiments that show plots of f+(t). They are [36, 37]. 
We will dicuss the data from [37] below. 

6.5 Ao and A + 

Most experiments have analyzed their data assuming linear form-factors. In Table 2 we 
have quoted the PDG2002 values and the more recent experiments not included in it. We 
will not use these numbers much, given the possible shifts when introducing a curvature 
in the analysis. 

7 Inputs 

As relevant combinations we have obtained in our earlier work [14] experimental values for 
Ry 2 leading to 

R%* = (0.22 ± 0.02) • 10" 3 (7.1) 

and 

L r 9 = (5.93 ± 0.43) • 10~ 3 (7.2) 
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Figure 5: The data on /+(i)// + (0) from K e ° 3 . The CPLEAR data are from Ref. [34]. 
The older data are for comparison, Buchanan 75 is [31] and Birulev81 is [30]. The latter 
reference has two distinct data sets. 



Process 


Ref. 


A + 


Ao 


K% 


[2] 


0.033 ±0.010 


0.004 ±0.009 


Ktz 


[%e 


0.0282 ±0.0027 


0.013 ±0.005 


K% 


[2] 


0.033 ±0.005 


0.027 ±0.006 


K°a 


[2]//e 


0.0300 ± 0.0020 


0.030 ±0.005 


K° e3 


[2] 


0.0291 ±0.0018 




Kts 


[2] 


0.0278 ±0.0019 




Kts 


[36] 


0.0293 ±0.0025 




K% 


[33] 


0.0321 ±0.0045 


0.0209 ± 0.0045 


K+z 


[37] 


0.0278 ± 0.0023 




K% 


[38] 


0.02748 ± 0.00084 





Table 2: The PDG averages for A + and Ao and the values from the most recent experiments. 
\xe means that lepton universality has been used in the measurement. The result of [37] 
is an update of [39] which was included in the PDG averages. The last result [38] is 
preliminary. 
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which used the estimate R V1 = —0.49 • 1CT 5 GeV 2 . The two resonance estimates for the 
Ry-y quantities done in [14, 40] via naive vector-meson dominance and the chiral inspired 
variety lead to basically the same estimate 

R%* « -4 • 10~ 5 GeV 2 . (7.3) 

For inputs for the other parameters we use our fits that include the latest data [41]. 
These are the p 4 fit, and fit 10 to 13 in Table 2 in [13]. This is a reasonable variation of 
the various input parameters. 

We use the PDG2002 mass values for all the particles involved and 

F n = 92.4 MeV , F K /F n = 1.22 . (7.4) 

The amplitudes for the decays are calculated with the mass of the charged kaon and 
the neutral pion. Those for K® 3 with the mass of the neutral kaon and the charged pion. 

The scale we use in all the coupling constants and the loop integrals is fi = m p = 
700 MeV. Almost all conclusions are done from experimental determinations of the various 
parameters so the results are /^-independent. 

8 Numerical Results 

8.1 Size of the pure loop contributions 

In Fig. 6 we show the results from the pure loop diagrams for f+(t). The different lines 
are for different choices of pion, kaon and eta masses. They give some indication of the 
size of quark-mass isospin breaking, but it does not include the enhanced effect discussed 
in [4, 5]. In Fig. 7 we show the equivalent results for f-(t). 

8.2 Comparison with Refs. [10, 21] 

At this point we should also compare with the calculation of [10]. In that paper numerical 
results are quoted in Eqs. (91-94). We agree, if we use input L[ and masses the same as 
theirs, well with their numerical expressions for A^ 6 op /_ (their Eqs. (92) and (94)) and 
with their numerical value for Ap 6 ° p /+(0) (their Eq. (93)). We do not agree with their 
expression for the t-dependence of A 1 ^ f + , our numerical values differ by roughly a factor 
of two. Given the good agreement with the other values this is rather puzzling. 

We have attempted to trace the possible sources of discrepancy, and the present con- 
clusions are: i) recalculating with all L\ = or only L r 9 = does not lead to agreement. If 
we use as input instead L\ = 0.0082 we get good agreement with f+(t), but it spoils the 
agreement for f-(t). ii) We have used a different split for the reducible and irreducible 
parts of the integrals and a different subtraction scheme. Comparing direct subresults is 
therefore rather difficult. We have also compared our results for the pion electromagnetic 
form-factor [14] with those in [21]. We have a small discrepancy for the real part there but 
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Figure 6: The pure loop diagram contributions to f+(t). Shown are the p 4 and the p { 
results for several different choices of the pion, kaon and eta masses. 
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Figure 7: The pure loop diagram contributions to f-{t). Shown are the p 4 and the p' 
results for several different choices of the pion, kaon and eta masses. 
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a rather large discrepancy in the imaginary part. The imaginary part as calculated in [14] 
would give a phase of a few degrees in [14] while those of [21] give a phase significantly 
above 10 degrees towards 500 MeV. The latter is not compatible with the expected per- 
turbative buildup of the phase from ChPT. Hi) The main expected source of discrepancy 
is the slightly different renormalization scheme. The quantity we disagree most strongly 
on has a strong cancellation for the value of the L\ used between the pure loop part and 
the L\ dependent part, possibly amplifying differences in the result. 

8.3 f+(t) and Comparison with the CPLEAR and KEP-PS E246 
data 

The numerical expression for the p 6 contribution with the C\ = and the L\ from fit 10 is 

TV 

= 0.01462 + 0.0896353t + 0.0006313t 2 + 0.3414t 3 [K%] 

= 0.01424 + 0.0840569t + 0.0071463t 2 + 0.3493t 3 [K&] (8.1) 

with t expressed in GeV 2 and it is valid in the range < t < 0.13 GeV 2 . 

We now compare our ChPT expression at order p 6 to the CPLEAR data [34]. The 
latter data are normalized to one assuming a linear dependence. It is therefore that the 
polynomial fits done in Sect. 6.3 added a normalization factor as well. We now perform 
a fit using the inputs for the L\ from fit 10. The other choices of the L\ give essentially 
similar results. So we fit the CPLEAR data to 

U(t) = ^- ) {l + ff\t) + ff(t)). (8.2) 

The effect of R+q of Eq. (4.6) goes into a + while -R+J gives the C\ part of c + and 
gives the C[ part of A + (c+ and A + are defined in Eq. (6.1)). 

Notice that the fitted value, using the input from the pion electromagnetic form-factor 
for R+2i gives a value for R+* in good agreement with the naive expectation. Notice also 
that the presence of a curvature does change the fitted value of the normalization by a 
little less than one %. A rather important change due to the inclusion of the curvature 
is the effect on the value of the slope. Notice that, using the ChPT expression and the 
curvature as determined from the electromagnetic form-factor leads to 

A + = 0.0170 ±0.0015. (8.3) 

This value comes from the ChPT in the following way 

A+ = 0.0283(p 4 ) + 0.0011(loops p 6 ) - 0.0124(C[) . (8.4) 

The p 6 correction is about 30%. The difference with the conclusions on A + of [10] is to a 
large extent due to their fixing the normalization at one. 
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Figure 8: The CPLEAR data together with the normalized ChPT result without the C, 
the last fit reported in Table 1 and the linear fit done by the CPLEAR collaboration. 



The CPLEAR data together with the normalized ChPT result without the G\ contribu- 
tion, the last fit reported in Table 1 and the linear fit done by the CPLEAR collaboration, 
is shown in Fig. 8. 

We have performed a similar exercise for the data from the KEP-PS E246 experi- 
ment. This is shown in Fig. 9. Again, both the linear fit and the one with the curvature 
fixed from the pion electromagnetic form-factor have a similar y 2 . The difference in nor- 
malization, relevant for \V US \ is about 0.6% and 

Rfl = -2.5 • KT 5 . (8.5) 

The value for is somewhat different from the value determined from the CPLEAR 
data, but still compatible with the resonance estimate. The value for the slope is 

A+ = 0.0214 ±0.0018. (8.6) 

Notice that in both experiments, KEP-PS E-246 and CPLEAR, we have neglected the 
experimental systematic errors. A possible discrepancy can only be put in after a full 
experimental analysis is performed. But again, both the normalization and slope are 
changed significantly from the linear fit case. 



8.4 The scalar form- factor fo(i) 

The scalar form-factor as we have shown above is important since it can be used to deter- 
mine the p 6 constants needed to evaluate /+(0). In this section we show numerical results 
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Figure 9: The KEK-PS E246 data together with the ChPT result, with fitted and 
the linear fit of the KEK-PS collaboration. 

for f (t) and A(t). We have used here the value of Fx/F n = 1.22. 

In Fig. 10 we show the function f (t) — /o(0) at order p 4 and for the various sets of the 
L\ of Ref. [13]. The value of /o(0) = /+(0) is discussed in the next subsection. 

As can be seen the convergence from p A to p 6 is quite good. These are the curves 
labeled "fit 10 p 4 " and "fit 10". Notice that all fits of the L\ done at order p 6 (fit 10-13) 
give basically identical results. The fit of the L\ at p A (labeled "p 4 fit") deviates somewhat 
but this we consider an artefact as discussed above. For comparison we have shown the 
part due to t(l — F K /F n )/(m 2 K — m 2 ). 

A good fit over the entire phase space < t < 0.13 (t in GeV 2 ) is given by 



The error from the values of the different sets of L\ is about 0.0013 at t = 0.13 GeV 2 . 

We have not attempted to do a fit to any of the data given the experimental situation 
on Ao- We would however like to point out that the predicted curvature in fo{t) is small 
but of the same order as in f+(t). As we saw above, this curvature made a rather large 
change in the measured value of A + . A similar effect in Ao can thus not be excluded and 
should be studied experimentally. 



A(f) 
A(t) 



0.25763t + 0.833045t 2 + 1.25252t 3 [K° e3 ] 
0.260444t + 0.846124t 2 + 1.33025t 3 [K+] 



(8.7) 
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L[ dependence 

0.12 | , 1 1 




0.02 0.04 0.06 0.08 0.1 0.12 

t [GeV 2 ] 

Figure 10: The form-factor f (t) — /o(0) with the C\ = 0. Shown are the cases for the 
neutral kaon decays for various sets of the L\ together with old current algebra result. 

8.5 The value of /+(0) and A(0) 

The results for /+(0) with C\ = (which is equivalent to A(0)) are shown in Table 3. 
The isospin breaking shown is only an estimate, we have calculated the K^~ 3 case with the 
masses mjc+ and m^o and K® 3 with m^o and m n +. Further work on including isospin 
breaking fully to two-loop order is in progress [42]. As discussed above, we consider the 
results with the L\ determined at p A order rather extreme. We have also investigated how 
A(0) varies if we vary the L\ according to the errors and correlations determined in [13]. 
For fit 10, the 68% CL error gives 0.00124 and for fit 11 it gives 0.00273. Notice that this 
latter set allows for a very large variation of L\. We take the latter 0.00273 as a sign of the 





^£3 




4 

p 


-0.02266 


-0.02275 


p 6 loops only 


0.01130 


0.01104 


p 6 -Li fit 10 


0.00332 


0.00320 


p 6 -Li fit 11 


0.00375 


0.00355 


p 6 -Li fit 12 


0.00216 


0.00189 


p 6 -L t fit 13 


0.00539 


0.00526 


p 6 -Li p A fit 


0.00891 


0.00863 



Table 3: The various contributions to A(0) = /+(0)| cr=0 — 1. 
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variation with the L\, notice that includes all the p 6 fits given above. As a conservative 
estimate of this error we take half of the p 6 loop contribution as error and add to it the 
error from the L\. We thus obtain 

A(0) = -0.0080 ± 0.0057[loops] ± 0.0028[L[] . (8.8) 

The value of /+(0) is related to A(0) via 

/ + (0) = 1 + A(0) - A (C[ 2 + C r u ) [m\ - ml)' . (8.9) 

7T 

A naive estimate of C[ 2 can be made from scalar dominance of the pion scalar form-factor 
(SMD) leading to 

C[ 2 | 5MD = -^|«-1.0 -KT 5 . (8.10) 
The other combination can in principle be estimated from Ao via 

A » = (f " l ) + < 2C « + c «» K + < + -4|sw • <*■") 

As an example take A = 0.009 ± 0.010 leading to 

2C[ 2 + C 3 r 4 = (-1.0 ±1.7) -10^ 5 . (8.12) 
Putting both estimates together gives 

/+(0)U « 0.0 ±0.1. (8.13) 

i 

Essentially a 1% precision on / + (0) requires a measurement of A to 0.001 (about 5%), 
assuming we can determine C[ 2 with the relevant precision from other sources [26]. 

The estimate of the p & corrections given in [4] is for the analytic contribution (propor- 
tional to (m s — m) 2 ) and contributes to the C\ dependent part 

/+(0)U w -0.016 ±0.008. (8.14) 
If we use this value we get for /+(0) the present best estimate 

/ + (0) = 1 - 0.02266 [p 4 ] + 0.01130 [p 6 pure loops] + 0.00332 [p 6 LJ] - 0.016 [p 6 [4]] 
±0.0057[loops] ± 0.0028 [LJ] ± 0.0080 [C[[4]] 
= 0.9760 ±0.0102, (8.15) 

where errors have been added in quadrature. This satisfies the bound /+(0) < 1 [43] (as 
quoted in [4]). Notice that this bound when combined with the other results here gives a 
constraint on the values of C[ 2 + Cg 4 and thus on a combination of the slope and curvature 
in the scalar form-factor. The total value and error on /+(0) awaits an experimental 
determination of the needed constants but Eq. (8.15) is our present best estimate. 
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The net order p 6 contribution is roughly one order of magnitude smaller, and with the 
same sign, as the p A one, due to the sizable cancellation between the calculation of the 
loop contributions and the estimate of the contribution from the p % constants. We have 
reevaluated \V US \ following the basic steps presented in Sec. 7.3 of [15] but using as input 
the previous expression (8.15) obtained for the K + decay case. The result should only be 
taken as preliminary since the full analysis including the isospin breaking terms at p 6 is 
still missing as well as an analysis of the effect of the curvature in the form-factor on the 
experimental value. It explicitly contains: the e.m. corrections at one-loop order and the 
next-to-next-to-leading order correction in the isospin limit. 

\V US \ = 0.2175 ±0.0029, (8.16) 

this value is indeed compatible with world-average value \V US \ = 0.2196 ± 0.0023 within 
errors. Notice that the theoretical errors are larger due to the uncertainty in the C\ and 
they are at the same footing as the actual experimental ones. Another recent evaluation 
of V us along the same lines can be found in Ref . [44] . 

9 Summary and Conclusions 

We have performed a calculation to two-loop order of decays in the isospin limit. As 
far as we have been able to check, this calculation agrees analytically with the earlier on 
in [10]. We agree with some of the numerical results of that work but not all. 

For K i3 decay measurements we have shown how the value of /+(0) needed for the 
determination of V us can be determined from the slope and curvature of f (t) which can be 
measured in K^. It is possible that additional information from the pion scalar form-factor 
near allows the measurement of the slope only to be sufficient [26] . 

We have presented a present best value for /+(0) based on an estimate of the p 6 con- 
stants C[ 2 and Cg 4 . It is clear that this can be further improved after the above measure- 
ments are performed. From the calculated best value of /+(0) we can deduce directly an 
improved estimated value of the CKM matrix element \V US \. 

As can be seen from our comparison with the data for f+(t) the presence of curvature 
can make a sizable impact both on the determination of the value of the form-factor at 
zero and on the slope. The large change in the value A + we found is entirely due to this 
effect and was compatible with estimates of the C\ involved. 
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